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ABSTRACT 

We present the full results of our 3-year long Submillimeter Array survey 
of protoplanetary disks in the Orion Nebula Cluster. We imaged 23 fields at 
880 /iin and 2 fields at 1330 //m, covering an area of ~ 6.5 arcmin 2 and containing 
67 disks. We detected 42 disks with fluxes between 6-135 mJy and at rms noise 
levels between 0.6 to 5.3 mJy beam" 1 . Thermal dust emission above any free- 
free component was measured in 40 of the 42 detections, and the inferred disk 
masses range from 0.003 — 0.07 M & . We find that disks located within 0.3 pc of 
9 l Ori C have a truncated mass distribution, while disks located beyond 0.3 pc 
have masses more comparable to those found in low-mass star forming regions. 
The disk mass distribution in Orion has a distance dependence, with a derived 
relationship max(M disfc ) = O.O46M (d/O.3pc) a33 for the maximum disk masses. 
We found evidence of grain growth in disk 197-427, the only disk detected at 
both 880 fim and 1330 fim with the SMA. Despite the rapid erosion of the outer 
parts of the Orion disks by photoevaporation, the potential for planet formation 
remains high in this massive star forming region, with w 18% of the surveyed 
disks having masses > 0.01 M & within 60 AU. 

Subject headings: circumstellar matter - - planetary systems: protoplanetary 
disks — solar system: formation — stars: pre-main sequence 
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Introduction 



Observations of the circumstellar disks that accompany young stars provide important 
constraints on the planet formation process. Their fundamental properties, like mass and 
size, play a significant role in the potential to form planets, and may determine the types 
of planets that can form. Our current understanding of disk properties has been expanded 



regions, Taurus-Auriga and pOphiuchus (e.g., 


Bcckwith et al. 


1990; 


Osterloh & Beckwith 


1995; 


Andre & Montmerle 


1994 


; lAndrews & Williams 


2005. 


2007a 


). The stars in these low- 



mass star forming regions evolve in relative isolation where only internal processes affect 
their disk properties and evolution, and their observations represent an important starting 
point to understanding single star and disk evolution. 

The vast majority of stars in the galaxy, however, were not created in regions like 
Taurus-Auriga or pOphiuchus, but in rich clusters that contain massive (> 8M ) O-type 
stars. Surveys of low-mass star formation show that 70-90% of stars within 2 kpc of the 
Sun formed in rich, embedded clusters, ~ 75% of which currently contain massive stars 



( jLada fc Ladall2003l ). Moreover, the mass distribution of molecular clouds is weighted toward 
the most massive clouds, which are more likely to form rich clusters containing O-type stars 



sec 



Evansl 119991 ) . There is even clear evidence that our Sun was born in a massive star 
forming environment. Studies of primitive meteorites have revealed they contain the decay 
products of 60 Fe, which can only be produced by a supernova, unambig u ously placing the 



Sun's formation n ear at least one massive star ( ITachibana et al.l 12006 
Gaidos et al.ll2009h . 



Krot et al. 2005 



In rich clusters, external influences on disk evolution become important and they can 



threaten the development and persistence of protoplanetary disks. The hig 



i den sity of stars 



can lead to enhanced probabilities of dynamical interactions (iBonnell et al.ll2003t). wh ile UV 
radiation from O-type stars can evaporate the disk material (1 Johnstone et al.lll998[ ). Such 
is the case in the Orion Nebula, which is the nearest, young, massive star forming region. 
Hundreds of youn g stars lie with the central 0.2 pc of the Orion Nebula Cluster (ONC) 
( [Hillenbrand! 1 1 9 9 71 ) . where the single, ~ 1 Myr old, 45 M Q , 06-type star, 9 l Ori C, dominates 
the radiation field of the region. Normally, protostars like those in Orion would be deeply 
embedded leaving them highly obscured throughout their formation. But the massive OB- 
stars of the Trapezium cluster, located at the heart of the Orion Nebula, have carved out 
a cavity with their radiation, clearing out the molecular cloud material along our line of 
sight, leaving a "blister HIP regi on. The fortuitous viewing geometry leaves the stars only 
slightly extincted (A v < 2.5 mag; iHillenbrandl 1 199 71 ) . allowing observational access into the 
molecular cloud at optical to infrared wavelengths. 
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The young stars in Orion were initially identified at optical and radio wavelengths as 



ioniz e d envelopes surro unding neutral condensations ( IL agues &: Vidall 1 19 79c iMoran et al. 



19821 ; iGaray et al.lll987l) and wer e later hypothesized to be circumstellar disks surrounding 
young stars by IChurchwell et al.l (119871 ) through centimeter wavelength VLA observations. 
Churchwell et al.l (119871 ) measured high mass- loss rates of M ~ 10~ 6 — 10 il^yr -1 , which 
they reasoned implied a large reservoir of neutral circumstellar material that would render 
the young stars invisible at optical and infrared wavelengths unless it was distributed in 
a disk-like geometry. Hubble Space Telescope (HST) images of the region spectacularly 
confirmed the disk hypothesis by imaging the externally illuminated objects and revealing 
they were circumstellar disks of dust and gas that ar e evaporating under the intense UV flu x 
from the most luminous star of the cluster, 9 l Ori C (lO'dell et al.lll993t lO'dell &: Wenlll994l ). 



The young objects were named "proplyds" by lO'dell et al.l ( 119931 ). an acryonym for 
PROtoPLanetarY DiskS, because they appeared remarkably different from isolated disk 
systems observed in low-mass star forming regions. To date, nearly 200 proplyds have been 
discovered in Orion through high resolution optical and near- I R observations w ith HST 
do'dell & Wcmglll99d : Bally et aI.lll998aL I2OO0I Ismith et al.lhoOol : iRicci et al.ll2008h . Most of 
the proplyds appear strongly ionized; they are surrounded by com etary shaped ionized gas 



cocoons, with brig ht heads facing 9 1 Ori C, and tails facing away ( IMcCullough et al.l Il995 



Bally et al.lll998al ). A small percentage of the disks are dark silhouettes , seen in extinction 
against the bright background nebula (see iMcCaughrean &: Q'delllll996l ). or embedded sil- 
houettes, projected against the glowing light from their own ionized cocoons. The pure dark 
silhouettes show no evidence of being photoevaporated and are likely not in the center of 
the Trapezium Cluster but only appear so in projection. 



Given the mass-loss rates of the disks, (M 



10" 



10 7 M yr 1 ), which are high 



enou gh to remove a MMSN (0.01 M ) in < 10 5 yr (IChurchwell et al.lll987l : iHennev fc O'Dell 
199JJ), it was questioned whether enough raw material could remain in the disks to allow 



planets to form. However, there is a critical radius surrounding each star, r g , within which 
the gravitational pote ntial of the star confines material to the sy stem for longer than planet 
formation timescales ( IJohnstone et al.lll998r lAdams et al.l 120041 ) . While it takes 0.1 to 1.0 
Myrs to remove the outer disk beyond r g ~ ^|rS where a is the sound speed in the gas, an 



(Adams et al. 


2004; 


Clarke 


2007) 



proplyds, then planet formation remains a possibility despite UV photoevaporation. 

Although extensive observations have been made of the Orion proplyds at optical, in- 
frared and radio wavelengths, none of these observations were able to yield disk masses. The 
disks are purely absorbing at optical and near-infrared wavelengths, and therefore, obser- 
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vations were only capa ble of placing lower limits on their column densities and masses, of 
5 x 1O~ 4 M = 0.5Mj up (IMcCaughrean et al.lll998l ). At radio wavelengths, the emission from 
the proplyds is dominated by free-free emission originating in the ionized gas cocoons, which 
swamps the weaker dust-disk emission. In most cases, disk masses can only be estimated 
from measurements of optically thin dust continuum emission at millimeter wavelengths. 
The millimeter spectral energy distribution (SED) of a disk originates from dust located in 
the cool, disk midplane. While the infrared observations are sensitive to only ~ 1% of the 
total dust mass (or dust volume), the millimeter observations are sensitive to dust located 
out to hundreds of AU, where the majority (>90%) of the mass lies. Low column densities 
at these radii make the emission optically thin, leading the disk luminosity to be propor- 
tional to the sum of emission from all dust grains. This yields a direct re lationship between 
sub millimeter flux density F u and disk mass (IBeckwith et al.l Il990l ) , providing one of 
the best ways for masses to be measured. Disk masses are often compared to a reference 
called the Minimum Mass Solar Nebula (MMSN=0.01 M Q ), which represents a lower limit 
on the amount of gas and du st that was present in the primordial Solar System disk when 
the planets began formation (jWeidenschillingl If977l ) . 



Since their discovery by HST imaging, the Orion proplyds have been surveyed by many 
millimeter wavelength interferometers in an atte mpt to measure the ir masses and infer 
their potential to form planets (BI MA at 3.5-mm. iMundv et al.l Il995l ; OVRO at 1.3-mm, 



mm; 



Bally et al.l 1998b; OVRO at 3-mm, Eisner fc Carpenter 



2006; CARMA and SMA at 1.3- 



Eisner et al.l 120081 ) but few disks were detected due to the limited sensitivity to dust 



emission at these long wavelengths. Even at millimeter wavelengths, the substantial levels 
of free-free emission swamps the thermal dust emission from the disks, making mass esti- 
mates difficult to impossible. The dust emission increases sharply with frequency, having a 
dependence of F u oc z/ 2 , and so the first successful detections of dust emission from the Orion 
proplyds came at CTfo-miUimeter wa velengths, 880 /iin observations with the Submil li meter 
Array (SMA) (jWilliams et al.ll2005l ). Four proplyds were detected by lWilliams et al.l (120051 ) 
having disk masses of 0.013-0.028 M , exceeding the MMSN, and showed the capacity for 
Orion disks to be truly protoplanetary. 



Building on the success of the pilot study observations by IWilliams et al.l (120051 ) we 
carried out an SMA survey at 880 /im and 1330 /im of protoplanetary disks in the Orion 
Nebula Cluster from late 2006 to early 2010. The goal of the survey was to measure the disk 
mass distribution in a rich cluster in order to the study the impact of massive stars on disk 
evolution. In this paper, we present the results of our complete survey, with observations 
taken at 880 /im and 1330 /im, covering 23 SMA fields and 67 proplyds. We targeted only 
HST-identified circumstellar disks, as there is no ambiguity about the geometry of their dust 
distributions. Two main results from the survey have been published to date, showing the 



-5 - 



inner 0.2 pc of the Orion Nebula lacks the largest an d most massive disks observed in low- 
mass star forming regi ons JMann fc Williams 2009a ). while the outer regions, beyond 1 pc 
from 9 1 Ori C, do not JMann fc Williams! l2009bl ~ Here, we present new SMA observations 
taken of disks located at intermediate distances (d=0.2-l pc) from 9 1 Ori C, to study the disk 
mass dependence on distance in the Orion Nebula Cluster. We describe the observations and 
data reduction methods in $2] and present the calculation of disk masses in We examine 
the dependence of disk mass on size and location within the cluster and also discuss the 
results and implications for planet formation in rich clusters in §H 



2. Observations 

Submillimeter interferometric observatio ns of 23 fields c ontaining 67 proplyds were con 



ducted with the Submillimeter Array (SMA; iHo et all 120041 ) between 2006 Dec 27 and 2010 



Jan 29 on Mauna Kea. We used the compact and extended array configurations of the SMA; 
see Tabled] for a journal of the observations. The compact configuration was chosen for ob- 
servations prior to 2009, to provide the best phase stability, and maintain sufficient resolution 
(~ 2.5" at 880 /mi) to distinguish individual proplyds. Data taken in 2009-2010 were ob- 
tained using the extended configuration of the SMA, which takes advantage of improvements 
to the array, particularly bandwidth doubling. 

The phase centers of the observations (Table [lj Figured]) were chosen to simultaneously 
maximize the number of proplyds imaged while minimizing contamination from the bright, 
nonuniform background (see Section 13.2 j) . Double sideband receivers were tuned to local 
oscillator (LO) frequencies of 340.175 GHz, 350.175 GHz or 224.170 GHz (882, 857 or 1330 
jttm, respectively: see Table 1). For observations taken prior to August 2009, each sideband 
provided 2 GHz of bandwidth, separated by ± 5 GHz from the LO frequency. Subsequent 
observations benefited from an increased 4 GHz of bandwidth per sideband, separated by 
± 6 GHz from the LO frequency. Each observation listed in Table [1] represents a full 8 hour 
track, except the last 6 fields (18-23), which were shared between pairs of fields because of 
the increased bandwidth. For the initial fields, we simultaneously observed the CO (3-2) 
and HCN(4-3) transitions by using an LO of 350 GHz. When it became clear that the 
line emission from the disk could not be distinguished from the extended molecular cloud 
emission, observations were switched to slightly lower frequencies (340 GHz) where both the 
atmosphere and receiver performance improves (see Tabled]). CO (3-2), HCN(4-3), CO (2-1) 
and 13 CO(2-l) were all detected but maps show they were indistinguishable from the cloud 
background and we do not discuss them further. 



Weather conditions for all observations were good, with < 2 mm precipitable water 
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vapor, or r(225 GHz) < 0.1 for 880 /mi observations and < 5 mm, or r(225 GHz) < 0.24 
for 1330 yum observations. The atmosphere was very stable during the nights and system 
temperatures ranged from 100-400 K. The total rms noise levels after calibration ranged 
from 0.6 to 6.0 mJy Table [T] summarizes the relevant observational information. 

The raw visibilities for each night were calibrated and edited using the MIR software 
package. Amplitude and phase calibration were performed through observations of the bright 
quasars J0423-013 and J0530+135. The compact array observations were interleaved between 
20-minute target observations of the proplyds, and 5-minute observations of each of the 
two quasars. We shortened the target integration to 15-minutes for the extended array 
observations to increase our monitoring of the quasars. Passband calibration was conducted 
with one of 3C454.3, 3C279, or 3C273. The flux scale was derived using measurements of 
Titan or Uranus as primary flux calibrators; we used the values derived by Mark Gurwell, 



available at http://smal.sma.hawaii.edu/callist/callist.html The flux scale is estimated to 



be accurate to ~ 15%. 

The calibrated visibilities were weighted by system temperature and inverted, then 
cleaned to generat e the synthesized continuum maps shown in Figures [2] & [3] using MIRIAD 



(ISault et al.lll995l ). All line transitions were edited out of the observations before producing 
the final images. The maps were created after eliminating wf-spacings shorter than 27 kA 
(physical baselines shorter than 23-m) to filter out uniform extended emission corresponding 
to size scales larger than 7.5". The size scale was chosen to preserve compact emission 
from the disks, while minimizing contamination from bright, extended cloud background. 
Simulations of the background (see Section I3~2]) confirmed that the 27 kA cutoff would resolve 
out the majority of extended emission and reduce the rms noise levels by factors of 2 to 3. 



Results 



The 23 SMA fiel d s included a total of 67 HST-identifie d dis ks from the catalo gs of 
O'dell fc Won~5 fjl996h . bally et al.l feoooh . ISmith et al.l J2005L and iRicci et al.l toO^ \ and 
are shown in Figure [TJ Continuum flux densities, F Q b s , were determined for each source by 
summing the emission within a 2a contour, and correcting for primary beam attenuation. 
The rms noise levels were determined using as much of the emission-free regions within the 
SMA primary beam. In total, we detected 42 of the 67 surveyed disks with signal-to-noise 
ratio's > 3. 



Before disk masses can be calculated, we need to separate the flux contributions from 
thermal dust emission from the disks, Fd ust , the free- free emission from the ionized cocoons, 
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Fff, and the background molecular cloud emission, Ft 



7 obs — -Pdust + F S + F h 



The following sections describe how Fg, Fb g , and Mdisk were determined. 



3.1. Free- Free Emission, F s 

The radio-submillimeter SEDs for the 42 disks detected at > 3a at 880 /im with the 
SMA are s hown in Figures H] and Q W e used VLA observations at centimeter wave- 



lengths from lGarav et all (119871 ) ; iFelli et al.l (119931 ) ; IZapata et all (120041 ) , and B IMA, OVRO, 



CARMA and SMA upper limits and detections at millimeter wave lengths (IMundy et al. 
19951 : iBallv et al.lll998bl : feisner fc Carpenter! lioo3 : feisner et al.lboosf ) to define the free-free 



spectrum and then extrapolated it into the submillimeter regime. Fits to the free- free emis- 
sion (F v oc i/ -01 ) and dust emission (F v oc v 2 ) are overlaid on the SEDs to show their 
relative contributions and contrasting spectral dependences. The long wavelength, 6 cm to 
1.3 mm, data show a flat spectral dependence consistent with optically thin emi ssion, z/ 1 



but with a range, highlig hted by the grey scale, which we attribute to variability (IFelli et al. 



19931 ; IZapata et al.l 120041 ) . We avoided observations taken at wavelengths longer than 6 cm 
(5 GHz) in this analysis, in order to avoid the turnover frequency, where the free-free emis- 
sion becomes optically thick and no longer follows a z/~ ai dependence. These SEDs show 
how substantial the free-free emission towards the Orion proplyds is, swamping the thermal 
dust emission from centimeter to millimeter wavelengths, and demonstrating why higher fre- 
quency, submillimeter observations are necessary to detect dust emission from these disks. 

The disks within fields 12 through 23 were not detected in previous surveys at radio 
to millimeter wavelengths and, with the exception of disk 072-135, none of the disks show 
signs of photoevaporation through their HST images, therefore their flux contribution from 
ionized gas emission is expected to be negligible at submillimeter wavelengths (F s <C 1 mJy 
). Ffj is listed in Table [2], and is the maximum level of free- free radiation extrapolated to the 
observing frequency. After subtracting off the free-free contribution to the proplyd fluxes, 
40 of the 42 detected proplyds had dust emission in excess of the ionized gas emission. The 
SMA emission detected towards two proplyds, 168-326NS and 180-331, is consistent with 
free-free emission and these are therefore listed in Table [3] instead of Table [2J 
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3.2. Background molecular cloud emission Fb g 



The Orion proplyds lie within a blister HII region which sits in front of a giant molecular 
cloud. The submillimeter emission from the clou d was mapped with the SC UBA camera on 
the James Clerk Maxwell Telescope (JCMT) by I Johnstone &: Bally see background 

of Figure [TJ The unfiltered background intensity is ~ 1-3 Jy per 15" SCUBA beam, which 
corresponds to a flux of ~ 10-30 mJy per compact array SMA beam, or a contribution of 
~ 5 — 15 x 10~ 3 M© to a typical disk mass. The SMA, however, resolves out much of this 
emission if it is extended. 

To characterize the effects of a strong background on the noise properties of each field, 
we performed simulations of the interferometric response to the extended cloud emission. 
The simulations were conducted using the array configurations and wf-tracks from the ob- 
servations and applying them to the corresponding positions in the SCUBA map. For each 
field, the SCUBA data were Fourier-transformed and sampled over the observational uv- 
tracks, inverted and cleaned. Spatially filtered maps of the background were made using 
Mf-spacings > 27 kA, to filter out the uniform extended emission on angular scales > 7.5", 
exactly as used to produce the final SMA maps described in £j2j The flux contribution from 
the background emission, Fb g , was determined by summing the emission within an SMA 
beam sized aperture placed at the position of the individual proplyd within the simulated 
maps. The values of Fh g are listed in Table [21 Usin g the clump mass spectrum in Orion 
(I Johnstone et al.ll200lL l2006t iJohnstone fc Ballyll2006l ). we estimate the probability of back- 
ground contamination from small clumps (7.5"-15") not detectable by SCUBA observations 
to be negligible, with « 10 -3 detectable clumps per 35" SMA primary beam. 



3.3. Disk Masses 

The submillimeter continuum emission is produced by dust grains in the disk, which 
absorb stellar UV and optical photons, and re-emit them at longer, infrared to millimeter, 
wavelengths. Submillimeter wavelength emission is sensitive to dust throughout the disk, 
where the majority (>90%) of the disk mass lies. Low column densities beyond ~20AU 
imply that the emission is mostly optically thin and the submillimeter fluxes depend directly 
on the total mass present in the disk: 

d 2 FJT) 

M disk = ^-4- 2 

The largest uncertainty in this calculation of disk masses lies in the value of the dust 
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mass opacity, k u , the ratio of e ffective dust cross sec tion to mass, which depends strongly on 
the dust grain characteristics f IBeckwith et al.lll990l ). 



We determined disk fluxes, F dust , according to Equation [TJ see Table [2J Disk masses 
were then calculated for the 40 Orion proplyds detected in significant dust emission using 
Equation [2j We used a distance to Orion of 400 pc, wh ich is based on recent measurements, 



inclu ding parallax observations of stars in the cluster, ( Menten et al. 



2007) and the orbital motion of binaries in the cluster ( Kraus et al. 



2007 



2007 



Sandstrom et al. 



20091 ). We used the 



Beckwith et al.l (119901 ) dust mass opacity, k v = 0.1(^/1000 GHz) cm 2 g~ 1 , which implicitly 
assumes a 100:1 gas to dust mass ratio. We u sed a dust temperature of 2 K, th e average 
for disks in Taurus-A uriga and p Ophiuchus ( [Andrews fc Williams! 12005k l2007al ) ; see also 
Williams et al.l (120051 ) . The final disk masses and upper limits for the non-detections are 
listed in Tables |5] and [3J 

The conversion of flux to disk mass is conventional and was chosen for the simplic- 
i ty of comparison with oth e r conti nuum studies, for example, in Taurus and p Ophiuchus 
([Andrews fc Williams! 120051 . l2007al ) . Longer wavelength observations of Orion disks with- 
out substantial free-free emission show grain growth to beyond millimeter sizes (Ricci et al. 
2010b, in press), as is also seen in Taurus and p Ophiuchus. The effect lowers the overall dust 
opacity, k v . Less direct evidence based on disk accretion rates and ages also suggest that 



disk masses may, in general, be un derestimated by between a factor of 4-8 (IHartmann et al. 



19981 ; lAndrews fc Williamsll2007al ). The disk masses derived here should be considered lower 
limits but this is somewhat mitigated by their comparison to a low estimate of the MMSN. 

The mass sensitivity of this survey depends on many factors, including the varying 
levels of free-free and background emission contributing to the observed fluxes, as well as 
the location of the proplyds within each field. We, therefore, derived our completeness level 
by performing Monte Carlo simulations to determine the fraction of sources that could be 
detected at > 3a as a function of mass, depending on the above-mentioned characteristics. 
The results reveal that our survey is 100% complete for disk masses > 0.0084 M & and 
50% complete for M d > 0.003 M (see also lMann fc Williamsl[2009a} ). 



3.4. 1330 /im observations 

Although most proplyds appear to be strongly ionized, there are a small number of 
silhouette disks which show no evidence of photoevaporation through their HST images and 
have negligible free-free emission levels at radio wavelengths. For these silhouette disks, we 
attempted to detect dust emission at a slightly longer wavelength of 1330 pm. The detection 
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of dust at more than one wavelength, on similar baselines, allows us to study the spectral 
behavior of the Orion disk emission and constrain their dust properties by comparing them 
with disks in Taurus and pOphiuchus. 

In the Rayleigh- Jeans (R-J) limit [hv <C kT), the Planck blackbody function has the 
dependence B oc z/ 2 , resulting in the submillimeter flux emission behaving as a simple power 
law in frequency, F v oc u 2+l3 . The power-law index, (3, is a function of the dust grain size, 
shape and composition. If a uniform shape and composition are assumed, knowledge of 
f3 can constrain the dust size distribution and expose whether the growth of particles has 
begun in the disk. The determination of /3 can be complicated if there is significant deviation 
from the R-J criterion or if there is a contribution of optically thick emission to the disk 
flux (IBeckwith et al.lll990l ). For these reasons, we empirically describe the submillimeter 
continuum emission as F v oc u a , with a related to /3, but with a minor corr ection (typically 
A/3 ~a few tenths), for optica lly thick emission from the inner disk (see IBeckwith et al. 



19901 : IBeckwith fc Sargentlll99lh . 



We observed two fields at 1330 pm, which were centered on silhouette disks 182-413 
and 114-426 (Fields 2 and 3), neither of which were detected (see Figure [7] for observations 
of disk 114-426). The proplyd 197-427, which lies in the primary beam of Field 2, was 
detected at 1330 jum, with a primary beam corrected flux of 17.4 mJy (see Figure [61 Table 
HJ). Disk 197-427 is the only proplyd in Orion detected at both 880 /im and 1330 jum with our 
SMA observations. The slope of its submillimeter emission, determined from a ratio of its 
1330 /im, 880 /im fluxes, is a = 2.8 ± 0.1. This slope is intermediate b etween that found for 
a typical Class II disk in Taurus- Auriga or pOphiuc hus with a ~ 2 ([Andrews fc Williams 



20051 . l2007al ; iRodmann et al.ll2006l ; iRicci et al.ll2010l). an d the ISM, which contains sub-pm 
sized dust particles and has an a ~ 4 (jPollack et al.lll994l ). In both Taurus and p Ophiuchus, 



and recently, in Orion (see Ricci et al. 2010b, in press), the systematic change observed in 
the submillimeter slope with disk evolutionary state is most readily interpreted as due to 
grain growth in the disks to at least millimeter sizes. Therefore, the slope of Orion disk 197- 
427 suggests dust grain growth may be underway in this disk. Its spectral slope is actually 
more consistent with a Class I Taurus disk than a Class II disk, which we speculate could be 
due to photoevaporation of the surrounding circumstellar envelope, leading to the premature 
emergence of a Class II disk in Orion. 

Disks 182-413, 174-414 and 183-419 were detected at 880 /im, but not at 1330 /im, 
so we can place constraints on their spectral slopes using their 3a flux upper limits at 
1330 /im. Fluxes, upper limits and a values for these disks are listed in Table HI The esti- 
mated spectral slopes for Orion di sks lie within with the range found for disks in Taurus or 



pOphiuchus, which show a=l-4 ([Andrews &: Williams! 120051 . l2007al ; IRodmann et al.l 12006 



Ricci et al 



20101 ). 



3.5. Near-Infrared Excess Stars 

Near-infrared emission in excess of stellar photospheric emission is often interpreted 



to be due to the presence of circumstellar disks around young stars. IVicente fc Alves 



(120051 ) noted that HST-observations were only able to identify disks around ~ 50% of the 
IR-excess stars in the Orion Nebula Cluster and suggested that the remaining disks were 
too small, <0.15" ^60AU in extent, to be resolved by HST imaging. In addition to the 
HST-identified proplyds which were the primary targets of our survey, there were also 43 
near- IR sources lying within the SM A fields. With the exception of one source, ID=253 



from iHillenbrand fc Carpenter! (120001 ). which was detected at ~3<r in Field 1, only IR-excess 
stars with HST-identified counterparts were detected by the SMA. The vast majority of the 
infrared-only disks were not detected by our observations, implying they are not only small 



in size, but also low in mass. 



4. Discussion 

4.1. Non-Detection of the Giant Silhouette Disk 114-426 

Surprisingly, the most prominent Orion disk, 114-426, an 1100 AU silhouette disk seen 
nearly edge-on in HST images, was not detected by the SMA at either 880 /im or 1330 /zm 
(see Figure [7]). The non-detection places a 3a upper mass limit of 0.012 M & on the disk 
at 880 fim, which includes a background flux correction of Fb g = -12.8 mJy (s ee Table EH) . 



The p revious disk mass upper limit based on a 1330 /im OVRO non-detection ([Bally et al. 



1998bl ) is < 0.033 M , which we re-calculated using the improved distance to Orion of 400 
pc, and the s ame T, k v as our study. By studying the extinction of disk 114-426 at multiple 
wavelengths, iMcCaughrean et al.l (119981 ) derived a lower limit on its mass of 5 x 1O~ 4 M . 



Together, these limits impose a disk mass on 114-426 between 5 x 10 4 -1.2 x 10 2 M (0.5- 
12M Jup ). 

The stringent mass upper limit of 1O _2 M on 114-426 is puzzlingly low considering its 
large extent. But the non-detection may not be surprising if we consider the background 
emission toward this source. Our SMA simulations, described in Section 13.21 demonstrated 
that the ridge of bright background emission that is clearly visible in Figure (U, dominates the 
noise in the sensitive interferometric observations of the region. Although we deliberately 
planned many of our observations to avoid the bright ridge, we made an exception for 114- 
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426. The high contribution from the background emission toward 114-426 is evident in the 
rms noise levels of both 880 /iin and 1330 /zm observations, which are higher than for pointings 
not located near the bright ridge (see Figured Tabled]). The estimated background flux 
within an SMA beam-sized aperture toward 114-426 is Fb g = -12.8 mJy, the largest absolute 
contribution toward any disk in our entire sample (see Tables [2j [3]). Therefore, the substantial 
noise levels towards disk 114-426 are not likely due to instrumental issues, but due to small 
scale inhomogeneities in the background emission towards this source. We, therefore, believe 
that disk 114-426 was not detected by our SMA observations because of its location near the 
strong, contaminating, background molecular cloud emission in the Orion Nebula. 

The lack of millimeter emission detected toward disk 114-426 has also been attributed 
to the g rowth of dust particles to large sizes th at no longer emit efficiently at these wave- 
lengths ( IB ally et al.l Il998at iThroop et al.l 120011 ). It is well-established that grain growth 



is a rapid process in protop lanetary disks (e.g., iD'Alessio et al.l l2001t IWilner et al.l 12005 



Dullemond fc Dominikl 120051 ) . However, grain growth must be balanced by collisional frag- 
mentation in order to m atch observations that reveal disks remain rich in small dust grains 
for up to several Myrs (IDullemond fc Dominikl 120051 ; iBirnstiel et al.l 120091 ) . It is also diffi- 
cult to reconcile the grain growth explanation with images of the disk taken at A=0 .4-4/xm 
flMcCaughrean et alill99~8l IThroop et allboOli Ishuping et aDl2003l : iRicci et al.lbo08h . which 
require a population of Atm-sized particles to generate. Studies of d isk extinction as a func- 
tion of wavelength (IMcCaughrean et al.l Il998l ; IShuping et al.l 120031 ) have even shown that 
most dust grains along the disk edges of 114-426 cannot be significantly larger than 4-5/im. 
Therefore, the primary reason for the non-detection of disk 114-426 by SMA observations is 
interference from the strong, background emission toward this source, not grain growth. 



4.2. Disk Masses and Sizes 

Theoretical models of disk photoevaporation predict that the most extended disks 
should experience the highest mass-loss rates because material at larger radii is more loosely 
bound to the embedded star and it provides a broader surface area for photoevaporation. 
Mass-loss occurs in the disks when the sound speed in the gas exceeds the escape velocity: 
c s = {kT I '/im^a > v esc = (GM+/R) 2, or for disk radii Rdisk ^ M+/T, where T is the temper- 
ature of the UV-heated disk gas and M* is the stellar mass. Escape velocities increase towards 
smaller radii, requiring the disk gas be heated to higher temperatures, T ~ 1/ Rdisk, before 
mass-loss can occur. Mass-loss rates drop two orders of magnitude, from M ~ 10~ 7 Mm yr" 1 



for 1 00 AU disks to ~ lO^MQyr" 1 for 20 AU disks f lJohnstone et al.l Il998l : lAdams et al. 



20041 ). removing the outer disks (>50-100AU) in 0.1-1.0 Myrs, while allowing the gravita- 
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tionally bound, inner (< 50 A U) disks to endure for >10Myrs flAdams et al.l 120041 ; I Clarke 
20071 ). The youth of 8 l Ori C (jHillenbrandl Il997l ) implies a photoionizing age < lMyr, and 
we, therefore, expect the outer disks of stars located near the Trapezium Cluster to be 
stripped away while the inner disks persist. 

To investigate the relationship between mass-loss and disk sizes, we plotted our calcu- 
lated disk masses versus radii for 31 sil houette disks with reso lved sizes from HST imaging 
in Figure El disk radii were taken fr om I Vicente fc Alyes ( 2005 ). We did not use the sizes of 
the unresolved Orion proplyds from Vicente fc Alvesl ( 2005 ) in this analysis because they do 
not have directly measured disk radii. The vertical dashed line in Figure [8] represents the 
resolution of the HST images, 0.15" ~ 60 AU at the distance of Orion. We expect photoe- 
vaporating disks to lose mass and grow smaller simultaneously, moving them from the upper 
right corner of Figure E] to the bottom left, an effect that is consistent with the observations 
shown. 

The diversity of initial disk properties within this sample obscures a clear-cut signature 
of the relationship between disk mass and size. For example, if we assume all the disks have 
similar surface density profiles, with E = S (r / -Ro) -1 , then the disk masses, M = 27rE -Ro-R, 
depend not only on disk size, R, but also on the surface density normalization, So- Disk 
masses and sizes are further complicated by their dependence on stellar mass, which sets the 
gravitational radius for how far photoevaporation can erode the disk on planet formation 
timescales, and on distance from 6 l Ori C, which also influences the mass-loss rates. Yet, 
despite the complexity of the relationship, we find evidence of a correlation between disk mass 
and size. The largest disks tend to be the most massive, and there is a lack of large, low-mass 
disks. The probability of correlation between disk mass and radius ranges from 95.4%-99.95% 
(2.0a — 3.5a; see Table EJ), which were calculated using various censored s tatistical tests that 
incorporate the 3<r upper limits for the non-detections (jlsobe et al.l 1 19861 ). 



Spatially resolved observations of many disks in Taurus- Auriga and pOphiuchus have 
revealed that the disk surface density profil e s quite uniformly follow a S ~ r _1 dependence 
within ~ 100-200 AU (Undrews et al.l 120091 : kndrews fc Williams! l2007bh . We also plot in 
Figured three diagonal dashed lines, which are comparisons to M^k for three normalizations 
of the surface densities, E Q , where £ = S (r/i? ) _1 . The uppermost normalization plotted is 
the surface density at 5 AU for a MMSN-disk, E 5 (MMSN 60 ), which we define here to contain 
0.01M© within R = 60AU: S 5 (MMSN 60 ) = 0.01 M /2vr(6O AU)(5 AU) = 47.2 gem" 2 . The 
remaining two profiles shown are a factor of \/To and 10 lower. The photoevaporation of 
disks should shift them to smaller sizes and masses along these surface density scalings. The 
normalizations plotted reveal that the majority of Orion silhouette disks have low surface 
densities when compared with the standard, MMSN. We also estimated the normalized disk 
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surface densities at 5AU, £5, for the Orion silhouette disks to compare them with Taurus 
and pOphiuchus disks (see Figure[9]). The histograms of log (£5), are s hown for Orion on top 
and Ta urus//) Ophiuchus below, with data taken for the latter set from lAndrews fc Williams 
(l2007bl ). The Orion distribution shows a broad peak near a surface density of ~ 20 gem -2 
at 5AU. The MMSN surface density at 5AU is shown by the dashed line in Figure [9j 
The distributions between the regions are similar, with many disks appearing to have lower 
surface densities than the standard of reference, MMSN. 



4.3. Dependence of Disk Mass on Distance from 9 1 Ori C 



All theoretical models of protoplanetary disk photoevaporation by O-stars agree in pre 
dieting that the most rapid disk erosion occurs for the largest disks located near 9 1 Ori C 



Jjohnstone et alill998l: Istorzer fc Hollenbach)ll999l : iRichling fc Yorkelboool : IScallv fc Clarke 



2001 



Matsuyama et al. 



2003 



Adams et al 



20041 ) . In Orion, the incident radiation field 



from the massive stars can be describe d by its ratio to the interstellar (Habing) field as 



Go = 13,000/c£ c (IParravano et al.ll2003l ). I Adams et al.1 (120041 ) calculated the gas tempera- 



tures of photoevaporating disks under various radiation fields (see their Figure 2) and showed 

1 

that the disk temperatures scale roughly with the radiation field: T ~ Gq. Therefore, at 
increasing distances from the massive stars, UV-heating of the disks declines, T ~ 1/d, 
allowing larger and more massive disks (see previous section) to survive photoevaporation: 



R, 



disk 



1/T ~ d. 



In order to examine the dependence of disk mass on location in the cluster, we plotted 
our calculated Orion disk masses versus their projected distances from 9 1 Ori C in Figure 
ITOl The most massive disks in the cluster are found at the largest distances from 9 1 Ori 
C, while the erosion of the upper end of the mass distribution towards smaller distances is 
clear from this figure. It is also interesting to note that all of the non-detections lie within 
0.3 pc from 9 1 Ori C. As approximately half of the survey sample were non-detections, we 
used censored statistical tests to calculate correlation probabilities. Results of the statistical 
tests support a correlation between Orion disk masses and projected distances from 9 1 Ori 
C, with high probabilities that range between 98.74% to 100%; 2.5-4.0 a (see Table |5]). We 
used projected distances because the true distances of the disks from 9 l Ori C are unknown. 
As a consequence, depending on their radial distance from 9 1 Ori C, the disks could in fact 
lie at larger distances from the central cluster, which would shift their positions in Figure 
[TOl to the right. However, such an adjustment is not likely to severely alter the trend of 
increasing maximum disk masses with distance. And in spite of the intrinsic diversity of 
initial disk properties, which should produce an observable spread in any real correlation, we 
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find compelling evidence of a dependence between disk mass and distance from an O-star in 
the cluster. 

The most massive disks are the most readily detectable in our SMA survey and are 
subject to the least observational bias. The maximum disk mass envelope in Orion is traced 
by the long, dashed line across the top of Figure [10J The trend of increasing maximum disk 
mass with respect to distance is obvious, indicative of a declining mass-loss rate. We derive 
a power-law relationship to fit the maximum mass envelope, which is described by: 



max(M dtsk ) = O.O46M (-^) a33 . (3) 

If we assume the disk mass distribution in Orion had no initial dependence on distance, 
we can estimate how long it would take to photoevaporate it to its present state. The 
maximum disk mass at d ~ 1 pc is 0.07 M©, and at d ^0.01 pc is 0.015 M©. Assuming a 
mass loss-rate of M ~ 10~ 7 M©yr -1 at d — 0.01 pc, it would take ~0.6 Myrs to truncate 
the upper end of the disk mass distribution from M disk = 0.07 M© to 0.015 M©, which is 
consistent with the photoevaporation lifetime of 8 1 Ori C of < lMyr. 

We empirically determined the distance, d, at which the MMSN disk populations within 
and beyond d differed with the greatest statistical significance. We found 8/14 disks at 
d >0.3 pc have masses > 0.01 M©, while only 11/53 disks within 0.3 pc have comparable 
masses. Fisher's exact test indicates a 96.2% probability (2a) that disks at large and small 
projected distances from 9 1 Ori C have different frequencies. This distance represents a sta- 
tistical boundary for circumstellar disk destruction by external UV photoevaporation in the 
Orion Nebula Cluster. Orion disks located at d < 0.3 pc have previously been foun d to have 



a mass distribution that is truncated at its high end (see lMann fc Williams! l2009al ). Beyond 
0.3 pc from 6 1 Ori C, heating by UV radiation is insufficient to drive appreciable mass-loss in 
the disks and we find the range in disk masses observed, M disk = 0.004-0.07 M©, is similar to 
the range exhibited by Class II disks in Taurus and p Ophiuchus. We, therefore, expect that 
Orion disks located beyond 0.3 pc from 9 l Ori C evolve similarly to disks in low-mass star 
forming regions, Taurus and p Ophiuchus, which do not experience UV-driven photoevapo- 
ration by nearby massive stars. Furthermore, it is likely that the initial distribution of disk 
masses in Orion was similar to that observed in Taurus- Auriga and p Ophiuchus, with the 
Orion disks subsequently photoevaporated from the outside inward, and surviving masses 
and sizes that are strongly dependent on their stellar masses and distances from 8 l Ori C. 
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4.4. Planet Formation Potential 

Despite the clear signs of outer disk destruction observed in Orion, the prospects for in- 
ner disk planet formation are not significantly different than in low-mass star forming regions. 
Overall in Orion, 31% (21/67) of the disks in our sample have masses greater than or equal to 
the MMSN (0.01 M Q ). For comparison, the fraction of Class II d isks with masses > 0.01 M P) 



in Tau rus and p Ophiuchus are ~ 37% and 29%, respectively ([Andrews fc Williams! 12005 



2007aj). The fraction of Orion disks with similar properties to the inferred initial conditions 
of our Solar System, with masses > 0.01 M Q within 60 AU, is 12/67 (~ 18%), slightly higher 
than the percentage found in Taurus-Auriga of « 13 % and comparable to the estimated 



fraction of stars with gas giant planets within 20 AU (jCumming et al.l 120081 ). Presumably, 
the lower mass, sub-MMSN disks can form planets less massive than Jupiter, and including 
these may raise the prospects for planet formation even further. While only a third of the 
total surveyed disks have masses comparable to a MMSN, all of the disks detected in our 
survey («57%) have masses >0.28 MMSN. 



5. Summary 

We have presented the results of our submillimeter interferometric survey of the 880 /im 
continuum emission from 67 young circumstellar disks in the Orion Nebula Cluster. These 
observations were taken to study the disk mass distribution in the region, and its dependence 
on distance from the most massive star of the cluster, 9 1 Ori C. We find that UV photoevap- 
oration has rapidly (<lMyr) eroded the outer parts of the disks located near this massive 
star, simultaneously reducing their sizes and masses. We present evidence that shows disk 
masses correlate with distance from 9 l Ori C. Circumstellar disk destruction by UV photoe- 
vaporation has a statistical limit at the distance of 0.3 pc from 9 1 Ori C; beyond 0.3 pc, the 
flux of UV radiation drops and may not adequately heat the disk gas to initiate mass-loss. 
Despite the hostile environment in the Orion Nebula, the disks appear very similar to those 
observed in sites of isolated, low-mass star formation, Taurus-Auriga and p Ophiuchus, in 
that we see: 

• Similar disk mass distributions after allowances are made for photoevaporation of the 
outer disk edges in Orion. 

• Grain growth in one Orion disk, 197-427. 

• A similar percentage of disks with M > 0.01M Q within R < 60 AU, which is also com- 
parable to the estimated fraction of radial velocity exoplanets. 
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We also observed, but did not detect, the largest known disk in Orion, 114-426, due 
to interference from the strong background molecular cloud emission toward this source. 
The overall potential to form planets in a rich cluster containing massive stars like Orion 
is comparable to that found in Taurus-Auriga and pOphiuchus, a promising result in the 
search for planets like our own. 
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Fig. 1. — Location of the Submillimeter Array (SMA) fields overlaid on a JCMT-SCUBA 
850 yum image of Orion. The black and white star near the center of the image marks the 
position of 6 l Ori C of the Trapezium Cluster and white crosses show the location of proplyds 
identified by HST observations. The solid white circles represent the 32" primary beam for 
SMA observations taken at 880 /im. Field numbers are labeled according to Table [TJ A 
dashed circle of radius 0.3 pc, or 155", around 6 l Ori C represents the sphere of influence 
of the massive star. Field 15, which contains proplyd 216-0939, is located a few arcminutes 
North of the image. 
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Fig. 2.— Aperture synthesis images of continuum emission at 880 /im taken with the 
Submillimeter Array towards proplyds the Orion Nebula Cluster. All proplyds detected at 
>3<t are shown and labeled with their names (see also Table [2]). The positions of the HST- 
identified proplyds are marked by crosses within the 10" x 10" field of view. Contours begin 
at the 3<r level, and each step represents 2a in intensity (see Table [2] for rms noise levels). 
The synthesized beam is shown in the bottom right corner of each map. 
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3. — Continued from Figure |2j 
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Fig. 4. — Spectral energy distributions for 42 proplyds detected at > 3cr with the Sub- 
millimeter Array (SMA) at 880 /im. The SMA measurements are represented by squares, 
millimeter observations by triangles and centimeter observations by circles. Open circles are 
upper limits from non-detections and uncertainties not shown are smaller than symbol sizes. 
The extrapolated range of optically thin free- free emission, F u oc u~ 0A , is overlaid in gray. 
A template to the disk emission, F v oc u 2 , is shown to guide the eye and reveal the relative 
contribution of the ionized gas and dust components. 
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Fig. 5. — Continued from Figure HI 
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Fig. 6. — Aperture synthesis images of continuum emission at 345 GHz (left) and 230 GHz 
(right) taken with the Submillimeter Array towards prominent silhouette disk 197-427 in 
the Orion Nebula Cluster. The rms noise, a, is specified in the bottom level corner of each 
map. The contours shown begin at the 3a level and each contour represents a step of 2a in 
intensity. The synthesized beam is also shown in the bottom right corner of each map. 
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Fig. 7. — Hubble Space Telescope (HST) and Submillimeter Array (SMA) images of disk 
114-426 in Orion. At left is an HST image of 114-426 taken from Massimo Robberto. The 
middle and right panels are aperture synthesis images of continuum emission towards disk 
114-426 at 345 GHz (middle) and 230 GHz (right) and show the disk was not detected at 
either wavelength. The rms noise, a, is specified in the lower left corner of the SMA maps. 
The contours shown begin at the 2a level and each step represents la in intensity. The 
synthesized beam is shown in the bottom right corner of the map. 



10 100 1000 

Disk Radius [AU] 



Fig. 8. — Disk masses versus radius for the 31 Orion silhouette disks with resolved sizes from 
Hubble Space Telescope (HST) imaging. Filled circles represent detections at > 3a, while 
open circles are the upper limits for the non-detections. The dashed vertical line represents 
the resolution limit of the HST images, 0.15" ^60AU. The diagonal dashed lines represent 
different normalizations of the disk surface density, which were calculated assuming E = S 5 
(r/5 AU)- 1 . The top normalization represents a disk mass of 0.01 M G (MMSN) within 60 AU, 
and the remaining profiles are a factor of vTO and 10 lower. 
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Fig. 9.— The distributions of the normalized disk surface densities at 5AU, log £5 in 
Orion (top) and Taurus/p Ophiuchus (bottom). For Orion, we show the surface densities for 
the sample of 31 silhouette disks with directly measured sizes from resolved Hubble Space 
Teles cope (HST) images. For th e Taurus/p Ophiuchus distribution, the data were taken 
from lAndrews &: Williams! ( j2007bl ). The dashed vertical lines represents the surface density 
at 5AU of the MMSN, which was calculated using S = S 5 (r/5AU) _1 , and shows the 
majority of disks in both regions have lower surface densities in comparison. 
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Fig. 10.— Circumstellar disk masses plotted against their projected distances from the 
massive star, 9 1 Ori C. Filled circles represent detections, while open circles are the 3a upper 
limits for the non-detections. The distance of 0.3 pc from 9 1 Ori C is marked by a dashed 
line, to separate the statistically different MMSN-populations. The maximum disk mass 
envelope is traced by the long dashed line across the top, to expose the absence of massive 
disks near 6 1 Ori C and the trend of increasing disk mass with distance. A derived power-law 
fit to the maximum mass envelope follows: max(M disk ) = O.O46M (d/0.3pc) a33 . 
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Table 1. Summary of Sub millimeter Array Observations 
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a 


(J2000) 


<5 (J2000) 
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UT Date 
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24 55.0 


340.175 


2008 Mar 14 


c 


0.09-0.11 


2.0 


2.6 x 1.5 


177.8 


9 (h) 


05 


35 21.70 


-05 


23 48.0 


340.175 


2008 Sep 17 


c 


0.05-0.06 


0.7 


1.8 x 1.5 


1.6 


10 (g) 


05 


35 17.70 


-05 


25 43.0 


340.175 


2008 Sep 24 


c 


0.05-0.15 


1.4 


2.0 x 1.6 


3.0 


11 (i) 


05 


35 17.10 


-05 


22 46.0 


340.175 


2008 Sep 30 


c 


0.05-0.15 


2.0 


1.8 x 1.5 


176.4 


12 


05 


35 29.40 


-05 


26 06.0 


340.175 


2008 Nov 11 


c 


0.06-0.08 


1.2 


2.2 x 1.8 


150.4 


13 


05 


35 07.20 


-05 


21 35.0 


340.175 


2008 Nov 26 


c 


0.08-1.0 


1.7 


2.1 x 1.6 


148.0 


14 


05 


35 13.20 


-05 


18 32.0 


340.175 


2008 Dec 21 


c 


0.06-0.07 


1.7 


1.9 x 1.5 


139.1 


15 


05 


35 21.60 


-05 


09 39.0 


340.175 


2008 Dec 23 


c 


0.10-0.15 


1.3 


2.1 x 1.6 


146.1 


16 


05 


35 25.30 


-05 


15 36.0 


340.175 


2008 Dec 24 


c 


0.13-0.16 


2.3 


1.9 x 1.4 


157.3 


16 


05 


35 25.30 


-05 


15 36.0 


340.175 


2009 Mar 26 


V 


0.05-0.10 


1.0 


0.3 x 0.2 


35.0 


16 


05 


35 25.30 


-05 


15 36.0 


340.175 


2009 Mar 27 


V 


0.03-0.05 


1.0 


0.3 x 0.2 


35.0 


17 


05 


35 18.10 


-05 


28 25.0 


340.175 


2009 Aug 25 


E 


0.04-0.06 


1.0 


0.9 x 0.7 


110.1 


18 


05 


35 34.67 


-05 


15 34.7 


340.765 


2010 Jan 19 


E 


0.02-0.03 


1.4 


0.7 x 0.6 


110.0 


19 


05 


35 28.04 


-05 


17 20.0 


340.765 


2010 Jan 19 


E 


0.02-0.03 


0.7 


0.8 x 0.6 


122.6 


20 


05 


35 05.40 


-05 


27 17.1 


342.000 


2010 Feb 28 


E 


0.06-0.08 


1.4 


0.8 x 0.7 


106.6 


21 


05 


35 32.10 


-05 


26 02.0 


342.000 


2010 Feb 28 


E 


0.06-0.08 


1.5 


0.7 x 0.5 


97.8 


22 


05 


35 12.10 


-05 


19 25.0 


340.765 


2010 Feb 29 


E 


0.05-0.08 


1.5 


0.7 x 0.6 


109.9 


23 


05 


35 14.05 


-05 


19 52.1 


340.765 


2010 Feb 29 


E 


0.05-0.08 


1.2 


0.8 x 0.7 


69.7 



Not e. — Notes — Column 1: Field Number; also labeled in Figure [T] Previous designation from l lMann fc Williams! 
2009b) listed in brackets. Column 2, 3: Phase Center Coordinates. Column 4: Observing Frequency. Column 5: UT Date 
of Observation. Column 6: Array configuration; C, compact (16-70 m baselines), E, extended (28-226 m baselines), V, very 
extended (68-509 m baselines). Column 7: Zenith optical depth at 225 GHz. Column 8: Root mean square noise measured in 
emission-free regions within the primary beam. Column 9: Dimensions of the naturally weighted synthesized beam. Column 
10: Position angle of synthesized beam, measured east of north. 
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Table 2. Disk Fluxes and Masses 



P roplyd 


r leia 


1 obs 


rms 


IT 

J? ff 


*bg 


r dust 


i\/r 

M disk 






[mjy] 


[mjy] 


[mjy] 


[mjy] 


[mjy] 


rn ni at 1 

Lu.uij.vi0 j 


w 




W 


(A \ 




(oj 


(<) 


K°) 


054-717 


20 


10.7 


1.4 


0.0 


0.0 


10.7 


0.52 ± 0.07 


072-135 


13 


18.6 


1.7 


0.0 


-2.8 


11.5 


1.05 ± 0.08 


121-1925 


22 


15.0 


1.5 


0.0 


0.0 


15.0 


0.74 ± 0.07 


132-1832 


14 


16.5 


1.7 


0.0 


0.2 


16.3 


0.80 ± 0.08 


136-1955 


23 


77.6 


1.2 


0.0 


0.0 


77.6 


3.80 ± 0.07 


141-1952 


23 


30.6 


1.2 


0.0 


0.0 


30.6 


1.50 ± 0.06 


155-338 


5 


39.6 


5.3 


4.0-14.0 


2.6 


23.0 


1.13 ± 0.34 


158-323 


7 


25.4 


2.4 


4.3-11.5 


-5.6 


19.5 


0.96 ± 0.21 


158-326/7 


7 


31.9 


2.4 


4.0-12.0 


-0.5 


20.3 


1.00 ± 0.19 


159-350 


5 


69.8 


5.3 


3.2-11.0 


0.3 


58.4 


2.88 ± 0.29 


161-328 


7 


38.7 


2.4 


0.1-1.5 


-1.2 


38.4 


1.89 ± 0.15 


163-249 


11 


34.8 


2.0 


0.0 


-1.6 


35.0 


1.72 ± 0.13 


163-323 


7 


22.1 


2.4 


2.4-3.0 


0.2 


18.9 


0.93 ± 0.15 


167-317 


7 


37.2 


2.4 


14.5-25.5 


-5.4 


17.0 


0.84 ± 0.18 


168-328 


7 


13.7 


2.4 


1.9-2.5 


1.0 


10.3 


0.50 ± 0.12 


170-249 


11 


14.6 


2.0 


0.0-2.5 


-0.6 


13.7 


0.68 ± 0.10 


170-337 


1 


19.1 


1.5 


4.0-9.0 


3.0 


7.1 


0.32 ± 0.08 


171-340 


1 


46.4 


1.5 


0.0 


-2.8 


49.2 


2.24 ± 0.08 


173-236 


1 


25.2 


2.0 


0.0 


2.9 


20.8 


1.03 ± 0.13 


174-414 


2 


12.7 


1.5 


0.0 


1.6 


11.1 


0.51 ± 0.09 


175-355 


8 


9.4 


1.5 


0.0 


-0.1 


9.5 


0.43 ± 0.11 


176-543 


10 


12.8 


1.4 


0.0 


-0.6 


13.4 


0.66 ± 0.07 


177-341 


1 


26.9 


1.5 


5.0-13.0 


-5.9 


19.8 


0.90 ± 0.07 


181-825 


17 


54.8 


1.0 


0.0 


0.0 


54.8 


2.70 ± 0.05 


182-413 


2 


13.3 


1.5 


0.0 


-3.0 


16.4 


0.74 ± 0.07 


183-419 


2 


6.2 


1.5 


0.0 


0.2 


6.1 


0.28 ± 0.07 


197-427 


6 


57.7 


1.9 


0.0 


-1.5 


59.3 


2.72 ± 0.10 


198-448 


8 


16.4 


2.0 


0.0 


-2.9 


19.3 


0.95 ± 0.12 


203-506/4 


8 


14.3 


2.0 


0.0 


-2.9 


17.2 


0.85 ± 0.13 


205-421 


6 


54.8 


1.9 


0.0 


-1.7 


56.5 


2.59 ± 0.10 


206-446 


8 


63.2 


2.0 


0.0 


0.0 


63.3 


3.12 ± 0.12 


216-0939 


15 


91.9 


1.3 


0.0 


0.0 


91.9 


4.53 ± 0.06 


218-339 


9 


6.6 


0.7 


0.0 


0.0 


6.6 


0.32 ± 0.04 


218-354 


9 


42.3 


0.7 


0.0 


-0.9 


48.1 


2.37 ± 0.04 


253-1536a 


16 


134.2 


1.0 


0.0 


-0.8 


135.0 


6.65 ± 0.05 


253-1536b 


16 


37.6 


1.0 


0.0 


-0.8 


38.4 


1.89 ± 0.05 


280-1720 


19 


21.8 


0.7 


0.0 


0.0 


21.8 


1.06 ± 0.03 


294-606 


12 


8.8 


1.2 


0.0 


-0.5 


9.3 


0.46 ± 0.06 


321-602 


21 


8.6 


1.5 


0.0 


0.0 


8.6 


0.42 ± 0.07 


347-1535 


18 


8.2 


1.4 


0.0 


0.0 


8.2 


0.40 ± 0.07 



Note. — Notes — Column 1: Proplyd designation based on the nomenclature of 
O'Dell & Wen (1994). Column 2: Observed Field in Figure [l] and Table [T] Column 
3: Integrated continuum flux density, corrected for SMA primary beam attenuation. 
Column 4: 1<t statistical error. Column 5: Range of extrapolated contribution of 
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free-free emission at 880 fim. Column 6: Estimated flux contribution from cloud 
background. Column 7: Derived dust continuum flux from the disk. Column 8: Disk 
mass (error does not include uncertainties in the flux scale of ~ 15%). 
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Table 3. Upper Limits on Disk Fluxes and Masses 



Proplyd 


r Ida 


3 a 


l ff 




M-disk 






[mjy] 


[mjy] 


[mjy] 


rn m A T 1 
[O.OIM0J 




U) 


(.•3 J 


(A \ 


I-'.) 


la\ 
K") 


114-426 


3 


< 10.5 





-12.8 


< 1.15 


159-338 


6 


< 17.4 


0.0-1.2 


-0.9 


< 0.85 


160-353 


6 


< 19.3 


0.0-2.0 


0.5 


< 0.82 


161-324 


7 


< 10.0 


2.5-4.5 


-0.8 


< 0.31 


163-317 


7 


< 12.8 


4.6-12.0 


-5.7 


< 0.32 


165-235 


11 


< 9.3 





1.6 


< 0.38 


165-254 


11 


< 8.0 





-3.7 


< 0.57 


166-250 


11 


< 7.1 





-0.3 


< 0.37 


166-316 


7 


< 11.6 


0.1-1.5 


1.7 


< 0.42 


167-231 


11 


< 10.5 





-0.6 


< 0.55 


168-235 


11 


< 8.4 





2.0 


< 0.32 


168-326NS 


7 


19.6 


8.6-18.0 


1.2 


< 0.02 


169-338 


1 


< 6.0 





0.7 


< 0.24 


169-549 


10 


< 4.6 





-1.3 


< 0.29 


171-334 


1 


< 5.7 


1.9-4.7 


1.2 


< 0.07 


173-341 


1 


< 4.6 





1.6 


< 0.14 


175-251 


11 


< 6.7 





-0.5 


< 0.36 


176-325 


7 


< 10.2 


2.8-4.0 


1.9 


< 0.21 


177-541 


10 


< 4.3 





-0.8 


< 0.25 


179-353 


1 


< 6.5 





1.3 


< 0.24 


179-534 


10 


< 4.7 





-0.2 


< 0.24 


180-331 


1 


7.7 


2.1-6.7 


2.1 


< 0.07 


181-247 


11 


< 9.3 


1.6-3.2 


0.9 


< 0.20 


182-332 


1 


< 6.1 





-2.1 


< 0.37 


183-405 


2 


< 5.9 





1.5 


< 0.20 


184-427 


2 


< 7.0 





1.7 


< 0.24 


213-346 


9 


< 2.6 





0.0 


< 0.13 



Note. — Notes — Column 1: Proplyd designation based on 
the nomenclature of O'Dell & Wen (1994). Column 2: Observed 
Field in Figure [l] and Tabic [l] Column 3: 3 a flux upper lim- 
its (or integrated continuum flux density for 168-326N, 180-331), 
corrected for primary beam attenuation. Column 4: Extrapo- 
lated contribution of free-free emission at 880^tm. Column 5: 
Estimated flux contribution from cloud background at 880/im. 
Column 6: 3 it disk mass upper limit. 
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Table 4. Submillimeter Continuum Slopes 



Proplyd 


Field 




Fl330fjm 


a 






[mJy] 


[mJy] 




(1) 


(2) 


(3) 


(4) 


(5) 


197-427 


2 


57.7 


17.4 


2.8 ± 0.1 


174-414 


2 


12.7 


<3.4 


> 3.0 ± 0.1 


182-413 


2 


13.3 


<3.0 


> 3.4 ± 0.1 


183-419 


2 


6.2 


<3.1 


> 1.6 ± 0.1 


114-426 


3 


<10.5 


<12.0 





Note. — Notes — Column 1: Proplyd designation 
based on the nomenclature of O'Dell & Wen (1994). Col- 
umn 2: Observed Field in Figure[l] Column 3: Integrated 
continuum flux density at 880 pra, corrected for primary 
beam attenuation. Column 4: Integrated continuum flux 
density at 1330 £tm, corrected for primary beam attenua- 
tion. Column 5: Submillimeter SED Slope as defined in 
Section 13.41 comparison of SMA and CARMA observa- 
tions. 



Table 5. Statistical Correlation Tests of Orion Disk Properties 



Correlaton Test 


Mass vs. Size 


Mass vs. Distance 




P(%) 


P(%) 


(1) 


(2) 


(3) 


Cox Hazard 


99.95 (3.5ct) 


99.99 (4.0a) 


Kendall's r 


95.36 (2.0ct) 


98.74 (2.5<r) 


Spearman's p 


97.17 (2.2<t) 


99.07 (2.6ct) 



Note. — Notes — Column 1: Censored statistical cor- 
relation tests were use d that incorporate 3cr upper limits 
for the non-detections l llsobe et al.lll986r ). These tests in- 
clude the Cox Proportional Hazard Model, the General- 
ized Kendall's t, and Spearman's p. Column 2: Probabil- 
ity of correlation between disk size and mass for silhouette 
disks. Column 3: Probability of correlation between disk 
mass and the projected distance of the disk from 8 1 Ori 
C. 
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